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The value of the tensor-to-scalar ratio r in the region allowed by the latest Planck 2015 measure¬ 
ments can be associated to a large variety of inflationary models. We discuss here the potential of 
future Cosmic Microwave Background cosmological observations in disentangling among the pos¬ 
sible theoretical scenarios allowed by our analyses of current Planck temperature and polarization 
data. Rather than focusing only on r, we focus as well on the running of the primordial power 
spectrum, as and the running thereof, /3s- R future cosmological measurements, as those from the 
COrE mission, confirm the current best-fit value for j3s / 10“^ as the preferred one, it will be 
possible to rule-out the most favoured inflationary models. 

PACS numbers: 98.70.Vc, 98.80.Cq, 98.80.Bp 


I. MOTIVATIONS 

The smoking-gun of inflation [1-3] is the detection of 
a stochastic background of gravitational waves. Such 
primordial signature is characterized by its amplitude, 
parametrized via the tensor-to-scalar ratio r. Recent 
analyses from Planck 2015 [4] have presented the tight¬ 
est bounds to date on r using temperature and polar¬ 
ization measurements. Albeit current Planck constraints 
are perfectly compatible with a vanishing tensor-to-scalar 
ratio, yet there is still enough room for other theoretical 
possibilities besides the Starobinsky i?^-gravity scenario, 
which emerges as the best-fit model. Looking forward to 
the next generation of CMB observations, and depend¬ 
ing on the value of r that Nature has chosen, one can 
envision two distinct possibilities: (a) either r turns out 
to be way too small to be measured by the next genera¬ 
tion of CMB observations, or (b) the value of r is large 
enough to be detected. However, in this latter case, the 
measured tensor-to-scalar ratio will typically correspond 
to several inflationary models. Given that measuring r 
(if r < few X 10“"^) might be extremely difficult [5, 6], 
and disentangling between the various models that lie in 
the same regions in the canonical (jis, r) plane might not 
be straightforward either, we explore here the possibil¬ 
ity of extending the analysis to other (complementary) 
inflationary observables. 

For the scalar power spectrum of the primordial per¬ 
turbations, we consider, as additional observables, the 
running and the running of the running j3s- For the 
primordial tensor power spectrum, we consider its run¬ 
ning nt- The aim of this paper is to assess the poten¬ 
tial of future CMB observations in falsifying inflation 
(or unraveling the fundamental model among the most 
favoured candidates after Planck 2015 data) by looking 
to these three additional observables. For illustration, we 
will consider some well-motivated models that are com¬ 
patible with current data. The structure of the paper 
is as follows. Section H deals with the basic definitions 


of the different cosmological observables and their cur¬ 
rent constraints. Section HI describes the theoretical 
predictions from the most favoured inflationary scenar¬ 
ios after Planck 2015 CMB temperature and polarization 
measurements. In Section IV we perform Markov Chain 
Monte Carlo (MCMC) analyses of the Planck 2015 data 
release. Our Fisher matrix forecasts in Section V show 
that, if the future preferred value of /3s is close to the 
current best-fit from Planck, future CMB probes may 
falsify the currently best inflationary scenarios. We shall 
conclude in Section VI. 


II. BASIC DEFINITIONS 


The power spectrum of the primordial curvature per¬ 
turbation, /, seeding structure formation in the universe 
is defined as 


iCitCp) = {27Tf S^^\p+k)Pdk) , ( 1 ) 


where the dimensionless amplitude of primordial pertur¬ 
bations ^({k) is defined through 

P^{k) = ^A^(fc). (2) 

The scale dependence of A^(fc) is parametrized by the 
spectral index: 


ris = 1 + 


dlnA2 

dln/c 


(3) 


Likewise, one can also define the scale dependence of the 
spectral index, which is called the running, as 


^ dus 
dlnfc ’ 


(4) 


as well as the running of the running, defined as 


Ps = 


dcts 

dlnfc 


(5) 
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In all these definitions, it is understood that quantities 
are evaluated at horizon exit fc* = aH {k^ = 0.05Mpc“^ 
throughout this study). In terms of the above parame¬ 
ters, the primordial power spectrum reads 

/ 7 \ rig —ln(/i:/fc+) + ^/3g \n^{k/k^) 

Alik) = Alik.) 

( 6 ) 

In the context of slow-roll, one can have a general idea 
about the magnitude of the above inflationary parame¬ 
ters in terms of the number of e-folds N. If we consider 
the empirical relation [7-9] — I oc 1/iV, one expects 

that 

“s ^ < 10""^ and /3s ~ ^ < 10"^, (7) 

for typical choices of the number of e-foldings = 50 — 
60. The latest Planck 2015 temperature and polarization 
TT,TE,EE-|-lowP [4] data analyses with r = 0 provide 
the following constraints: 

Us = 0.9586 ± 0.0056 , 

Qfs = 0.009 ±0.010 , 

/3s = 0.025 ±0.013 . 

What is interesting to notice in these constraints, is a 
slight preference for a positive /3s ^ 10“^, while as we 
will explain shortly, slow-roll inflation predicts typically 
a smaller and negative Ps. 

The tensor contribution to the primordial power spec¬ 
trum is parametrized by the tensor-to-scalar ratio r 

r = Ptik.) / P(ik.) , (8) 

where Ptik) = is the tensor power spectrum, 

and it is parametrized at first order as 

Alik) = Alik.) , (9) 

in which nj is the spectral index of tensor modes. In 
the slow-roll regime, the magnitude of r can vary within 
a large range, and this is the main difficulty in test¬ 
ing inflation through the detection of B-modes. This 
can be understood in the context of phenomenological 
parametrizations of inflation [7-9]. In such approaches, 
the (rig, r) plane appears to be unevenly filled, and one 
can even argue on the existence of a “forbidden zone” 
in the r-direction, depending on the precise value of rig, 
see Figure I. Future CMB missions aim to reach the im¬ 
portant theoretical milestone of r = 2x I0“^-(60/iV)^[I3], 
which would signal super-Planckian inflaton excursions 
[14-16]. 


^ This observation has been made previously in different contexts 
in [10-12]. 


III. MOST FAVOURED INFLATIONARY 
SCENARIOS 

In the following, we shall review the most favoured 
models (including their predictions for the different in¬ 
flationary observables: r, rij, rig, ag and /3g) after Planck 
2015 data release. 


A. Quadratic scenarios 

This class of scenarios represents the simplest theoret¬ 
ical possibility. It includes: 

The chaotic scenario, V oc (jP, both with mini¬ 
mal and non-minimal coupling to gravity [17-23]. The 
former is disfavoured with respect to the latter so the 
non-minimally coupled version is perfectly compatible 
with current data [24]. The predictions in the (ng,r), 
(ng,Q;g), (ng,/3g) and (nt,r) planes for these two models 
and ^TZcj)'^) are depicted in Figures 1, 2, 3 and 4 for 
two possible choices for the number of e-folds, iV = 50 
and N = 60 Notice, from Figure 1, that the trajecto¬ 
ries in the (ug, r) plane for the non-minimally coupled 
case i^TZpP) start always at the point corresponding to 
the (jP model predictions and then, as the coupling ^ 
takes positive values, the tensor contribution is reduced, 
and the scalar spectral index rig is pushed below scale 
invariance, see Ref. [24]. Negative values of the coupling 
^ (not illustrated here) are highly disfavoured by current 
CMB observations, since they will lead to large values 
of the tensor-to-scalar ratio r. Concerning the running 
of the scalar spectral index Og, the trajectories for the 
two quadratic scenarios considered here are depicted in 
Figure 2. Notice that positive values of the coupling ^ 
will change the predicted value of ckg in the (jP scenario 
(ofg = corresponding to ccg = —0.00056 for 

N = 60) to slightly larger values, albeit the trajectories 
always stay in the ag < 0 sub-plane. The running 
of the running parameter, /3g, barely changes with 
respect to its predicted value in the non-minimally 
coupled case (i.e. ^ = 0, for which /3g = —4/V^, giving 
/3s ~ —1.8 X 10“® for N = 60) as the coupling ^ gets 
positive values, see Figure 3. Finally, in Figure 4 we see 
that all models follow the theoretical curve nt = —r/8. 
In particular, the chaotic (j)^ model predicts a tensor 
spectral index of nt ~ —0.019 (uj ~ —0.016) for N = 50 
iN = 60); an increasing positive value of within the 
non-minimally coupled model, diminishes the predicted 


^ The value of ^ ranges from ^ = 0 to ^ = 0.0065 in Figures 1, 2 
and 3. 

® The case of ^ = 0 is equivalent to the standard inflationary 
chaotic scenario in which the predictions are n. = 1 — 2/N and 
r = S/N, corresponding to Us = 0.967 and r = 0.13, respectively, 
for N = 60. 
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value to nt — —0.018 {rit — —0.010) for ^ ~ 0.0059 and 
N = 50 {N = 60). 

The Natural inflation scenario (minimally cou¬ 
pled to gravity), where the inflaton is a Pseudo-Nambu- 
Goldstone-Boson (PNGB), which potential is invariant 
under the shift (p ^ (p + 27r/, and it is given by 


V{cP) = Po [1 - cos(<^//)] , (10) 


with / the PNGB decay constant [25-27]. It is straight¬ 
forward to perform the slow-roll analysis and obtain 
the analytical expressions of the spectral index and the 
tensor-to-scalar ratio: 


■e“^ (1 + a/2) -h 1' 


(1 + o/2) - 1 ’ 

where the parameter a is defined as ^ a = No¬ 

tice that for small a (i.e. very large values of f) the pre¬ 
dictions of the natural inflation scenario coincide with 
those of the minimally coupled chaotic inflation model 
V (X (p'^. Even if the flatness of the PNGB potential is 
protected by the shift symmetry, it is not clear whether 
this structure can be UV completed. For a recent discus¬ 
sion on the issue and some solutions see e.g. [28-32]. 

Figures 1, 2, 3 and 4 show the predicted trajectories in 
the {ns,r), (rig, Ug), (rig, /3g) and {nt,r) planes for N = 50 
and N = 60, and / varying from 3.45Mpi to lOOMpi. For 
the smallest value of / considered here, / = 3.45Mpi, 
a very small value of rig ~ 0.9152 is found. Agreement 
with Planck data implies that the decay constant satisfies 
/ > 5.3Mpi, for TV = 50 — 60. Larger values of / increase 
the value of the tensor-to-scalar ratio, until the prediction 
reaches the one of minimal chaotic inflation, as shown in 
Figure 1. In Figure 2, we illustrate that large values of 
/ lead to small values for the running of the spectral 
index, which eventually will reach the predictions for the 
minimal chaotic scenario. In contrast, the value of /3g, 
barely changes when / varies, remaining around in /3g ~ 
-3 X 10-5 and /3g ~ -1.7 x IQ-^ for TV = 50 and iV = 60, 
respectively, see Figure 3. Goncerning the tensor spectral 
index, for a value of / = lOOMpi the predictions coincide 
with those of the cp"^ model. Whereas lower value of /, 
corresponds to smaller values nt- For instance, nt ~ 
-0.0006 {nt ~ -0.0002) for / = 3.45Mpi and N = 50 
{N = 60), following the consistency relation n* = —r/8, 
as expected (see Figure 4). 


(1 + Q!/2) - 1 

8a 


( 11 ) 


B. Higgs-like scenarios 

This class of scenarios is described by a symmetry 
breaking potential, 



alike to the one of the standard model Higgs particle, 
but with a non-minimal-coupling to the Ricci scalar, 
see Refs. [33-35]. It also includes, as a limiting case (for 
^ —> oo), the i?^-gravity Starobinsky scenario [36]. Notice 
as well that the limiting case ^ 0 corresponds to the 

quartic potential scenario, V oc pP. One can find a suit¬ 
able set of inflaton potentials for different values of the 
inflaton vacuum expectation value v [23]. In this work we 
illustrate the predictions of a Higgs-like scenario for u = 0 
and for different positive values of as well as for N = 50 
and TV = 60 e-folds 5. Figure 1 clearly shows that the 
limiting case ^ > 0, corresponding to the quartic poten¬ 
tial is not in good agreement with Planck data, as its 
predictions for the inflationary parameters (rig ~ 0.941, 
r ~ 0.31 and Tig ~ 0.951, r ~ 0.26 for TV = 50 and 
TV = 60, respectively) are highly disfavoured. When the 
non-minimal coupling to gravity, is increased, the ten¬ 
sor contribution is reduced, while the predictions reach 
those corresponding to the Starobinsky scenario, as long 
as ^ > 10^. In this limit, rig ~ 0.961, r ~ 0.0041 
(rig ~ 0.968, r ~ 0.0023) for TV = 50 (TV = 60), val¬ 
ues which are in excellent agreement with current CMB 
data. 

Goncerning the running of the spectral index, increas¬ 
ing the value of ^ will drive the values of ag from the 
one corresponding to the quartic potential to slightly 
larger ones, corresponding to the Starobinsky scenario, 
keeping always the trajectory in the ag < 0 sub-plane 
(see Figure 2). The predictions of the running of the 
spectral index for the quartic (Starobinsky) scenarios 
are ag ~ -0.0011 (ag ~ -0.00074) for TV = 50, and 
ag ~ —0.0008 (ag ~ —0.00052) for TV = 60. As in the 
case of the previous models, the running of running of 
the spectral index, j3s, remains almost constant as ^ is 
varied, as shown in Figure 3. In particular, in the Higgs- 
like scenario, /3g ~ —3.5 x lO-^ (/3g ~ —2.5 x lO-^) for 
TV = 50 (TV = 60). This model allows for a wide range 
of values for the tensor spectral index, starting from 
the predictions from the model around nt — —0.039 
{nt — —0.039). Then, an increasing value of ^ pushes 
down the predictions for rij down to very small val¬ 
ues around nt — —0.0005 {nt ~ —0.0003) for TV = 50 
(TV = 60), (and thus coinciding with the values predicted 
from Starobinsky inflation), along the theoretical curve 
nt = —r/8 depicted in Figure 4. 


^ .^pl = I/VSttGjv — 2.43 x 10^® GeV is the reduced Planck mass. 


® The coupling A cancels out in the slow-roll calculations. 
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C. Hilltop scenarios 

For completeness, we should also consider this class of 
scenarios, described by potentials 

V{4>) = K) [1 - , (13) 

since its predictions in the (ns,r) plane lie very close to 
the ones associated to the models discussed before [15]. 
Within these scenarios, we can distinguish two sub-cases: 

1. p = 2, corresponding to the quadratic hilltop sce¬ 

nario, where inflation takes place close to a local 
maximum; = 0 and V"{(j)) < 0. 

2. p > 2, corresponding to a generalization of the sim¬ 

plest quadratic case, where here inflation happens 
close to a local maximum where additionally, higher 
derivatives of the potential vanish, i.e. = 

= • • • = = 0 and, again, 

< 0 . 

We restrict our analysis to the first case, p = 2, in 
which the spectral index and the tensor-to-scalar ratio 
read as 

ris = 1-4|77 o| 

with [ 770 1 = In Figure 1 we depict the predic¬ 

tions for this model in the plane (ns,r). The parameter 
770 varies from 770 = 10 “"^ to 779 — 2 x 10 “^, pushing 
to smaller values as 77 decreases. With 770 ~ 8 x 10“^ we 
obtain a tensor-to-scalar ratio of r = 0.0375 for the case 
N=60, and r ~ 0.0516 for N = 50, both corresponding 
to a spectral index ~ 0.968. Notice from Figures 2 
and 3 that for the same value of 770 ~ 8 x 10“^ we obtain 
a running of the spectral index of as — —0.00107 and a 
running of the running /Ig ~ —0.000065 (og ~ —0.00073 
and /3g ~ -0.0000386) for iV = 50 (TV = 60). In Figure 4 
we observe that this scenario predicts almost negligible 
values of the tensor spectral index for the range of val¬ 
ues of 770 commented above. The predictions reach the 
smallest values of rit found in this work: Ut — —0.0003 
{nt ~ -0.00006) for IV = 50 (IV = 60) and 770 ~ 2 x IQ-^. 

IV. CURRENT CONSTRAINTS 
A. Cosmological data and methodology 

We consider the new data on CMB temperature and 
polarization measured by the Planck satellite [37-39]. 
We use the Planck TT temperature-only likelihood (here¬ 
after Planck TT) and the Planck TT,TE, and EE power 
spectra data (hereafter Planck TTTEEE) up to a max¬ 
imum multipole number of £max = 2500 combined with 
the Planck low-£ multipole likelihood that extends from 
I! = 2 to £ = 29 (denoted as lowP). We use the Boltzmann 


Parameter 

Physical Meaning 

Prior 


Baryon density 

0.005 ^ 0.1 

uJc = 0,ch^ 

Cold dark matter density 

0.01 -> 0.99 

Os 

Angular scale of recombination 

0.5 ^ 10 

T 

Reionization optical depth 

0.01 ^ 0.8 

ln(10i°Ag) 

Primordial scalar amplitude 

2.7 -> 4 

ris 

Scalar spectral index 

0.9 ^ 1.1 

Oig 

Running of Tig 

-0.04 0.06 

Ps 

Running of as 

-0.04 -> 0.08 

r 

Tensor-to-scalar ratio 

0^2 


TABLE I. Uniform priors for the cosmological parameters 
considered in the present analysis. 


code CAMB [40] and generate MCMC chains using the 
publicly available package cosmomc [41]. We consider a 
ACDM extended model, described by the following set of 
parameters: 

{wbjWc, 0s,T,ln(lO^°Ag),77g,r, Q:g,/3g} . (15) 

In Table I, the uniform priors considered on the different 
cosmological parameters are specified. We do not con¬ 
sider the spectral index for tensor perturbations nt as 
an additional parameter in our MCMC analyses, since, 
as recently shown in [42], the current and future error 
bars on this parameter are considerably larger than the 
predictions of the different theoretical scenarios explored 
here. Therefore, the tensor spectral index is fixed in 
what follows to the slow-roll consistency relation value, 
nt = -r/8. 


B. Results 

While the latest Planck data provide evidence against 
some of the models explored here [4] , these measurements 
can not single out the responsible mechanism for the in¬ 
flationary process, nor to falsify this theoretical scenario 
by themselves. 

This can be noticed from the contours shown in Fig¬ 
ures 1 2 and 3, where it is clear that all the models 
described above have some trajectories in the (7ig,r), 
{ns, Us) and (7ig,/3g) planes which lie within the current 
68 % and/or 95% CL allowed regions. Figure 1 depicts 
the current 68% and 95% CL allowed contours in the 
(77 g,r) plane from Planck TT plus lowP data, as well 
as from Planck TT plus lowP data plus TTEETE mea¬ 
surements, together with the predictions from Natural, 
Hilltop, Higgs-like, quartic, chaotic ® and Starobinsky 


The chaotic model is studied both in its minimally and non- 
minimally coupled versions. 
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Parameter 

Planck TT-|-lowP 

Planck TT,TE,EE±lowP 

r (95% CL) 

< 0.27 

< 0.23 

Us 

0.959 ± 0.008 

0.9591 ± 0.0056 

0's 

0.0081 ± 0.014 

0.0077 ±0.011 

Ps 

0.034 ±0.016 

0.0313 ±0.014 


TABLE II. 95% CL constraints on the tensor-to-scalar-ratio r and mean values (together with their 68% CL errors) of ris, as 
and Ps obtained with the two possible data combinations considered in this study. 


inflationary scenarios, for both iV = 50 and = 60 e- 
folds. The addition of EE and TE spectra to Planck TT 
plus lowP data helps in constraining the scalar spectral 
index Us, however there is only a mild improvement in the 
tensor-to-scalar ratio upper bound. Notice, as previously 
stated, that the predictions for the inflationary parame¬ 
ters Us and r from these models are all well within the 
current 68% and/or 95% CL allowed regions and there¬ 
fore all of them (except for the case of the (/>'* potential 
with N = 50) are still feasible. One could ask if current 
measurements of other inflationary parameters, as the 
running of the scalar spectral index Og and/or its run¬ 
ning, /3s, may help in disentangling among the plethora of 
models still allowed by current data. Figure 2, illustrates, 
together with the trajectories in the (rig, Og) plane for the 
models explored here, the 68% and 95% CL allowed re¬ 
gions from Planck TT plus lowP data as well as from 
TTEETE plus lowP measurements. Notice that current 
bounds on as are unable to discard any of the possible 
inflationary models. Figure 3 shows the equivalent but 
in the (rig,/3s) plane. Interestingly, Planck measurements 
of /3g seem to exclude the value /3s = 0 at the ^ 2a level. 
The theoretical scenarios illustrated here could be ruled 
out with a much higher significance if the value of /3g 
preferred by Planck 2015 measurements (i.e. /3 ~ 0.025) 
is confirmed by future CMB data. We shall explore this 
possibility in the next section. 

Table II shows the 95% CL bounds on the tensor-to- 
scalar-ratio r as well as the mean values and 68% CL 
errors of the remaining inflationary parameters rig, as 
and /3g obtained with the two possible data combinations 
considered in this study. Notice that the limits on r are 
considerably relaxed when adding the running and the 
running of the running as additional parameters in the 
analyses. The mean values and the errors on rig and 
/3g are in very good agreement with those found by the 
Planck collaboration and reported in Ref. [4]. 


V. FORECASTS 


The aim of this section is to forecast the potential of 
future CMB satellites in constraining the {r, rig, Og,/3g} 
parameter space via the Fisher matrix formalism. 


A. CMB Likelihood 


Assuming that the fraction of sky surveyed /sky is the 
same for CMB temperature and polarization measure¬ 
ments, the likelihood associated to a single frequency 
CMB experiment can be written as 


- 2 In = ^(2^-b 1)/, 


sky 


( /^BB \ /yBB 

In ( ^ 1 - 


cr) 


c. 


BB 


C /^TT^EE f/^TE\2\ 

_WC L 1 _ 3 (16) 

CfTcEE _ ^cJEy j ^ ) 


aiTT/-iEE I /^TT/^EE eyy^TE/^TE 
j- 

(jTTfjEE _ (^(JTE'^2 


where the Cf-^ refer to the theoretical (mea¬ 

sured) power spectra for X,Y = T,E,B. Due to the 
finite resolution of the spectra, there will be an induced 
noise in the map that should be added to the Cg. In ad¬ 
dition, following [43] we will also include the foreground 
contribution to the map as a residual noise, and therefore 


Cg = Cf + Ng + Rf , (17) 


where will be our theoretical power spectra (com¬ 
puted by the Boltzmann solver codes CAMB [40] or 
CLASS [44]), Ng is the instrumental noise (which is a 
function of the frequency channel, see below) and Rf 
refers to the residual foreground subtraction (which will 
also depend on the frequency channel). This latter quan¬ 
tity reads as 


Rf{E) 


^ fore 


I 


Ng{n) 


4 Cliu) \ 

Nchan{Nchan — 1 ) CI{vf) J 


(18) 


where the first term corresponds to the uncertainty of a 
given foreground at a given frequency n, C}(iy) and ai(iy) 
represent the power spectra and the foreground subtrac¬ 
tion level, respectively. The second term in Eq. (18) takes 
into account for the instrumental noise of the channel 
at which the foreground model is constructed, and i/p 
is the frequency at which the foreground is modelled. 
In the case of a multifrequency experiment, as Planck 
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or COrE, the expression for the likelihood Eq. (16) still 
holds. However, in such a scenario, the total noise power 
that should be added to the Ci is written in terms of 
a weighted combination of the noises from the different 
channels [43]. Therefore, for a multifrequency experi¬ 
ment, Eq. (17) reads as 

Cl, = + Nf , (19) 

where the effective noise term is given by 

(7Vf)-2= ^ [{Rj{ui) + Ni{u,)) X (20) 

i,j>i 

{Rf (fo) + 2 ^ ■ 

We focus here on the future satellite experiment 
COrE [13], covering 70% of the sky. In the next sections 
we will describe the modelling of the experimental resolu¬ 
tion and the main foregrounds for this future CMB mis¬ 
sion, and therefore in what follows the numbers quoted 
will always refer to /sky = 0.7. 


mission ‘ and their up-to-date modelling, as provided by 
the Planck team. 


a. Synchrotron emission 


Synchrotron emission results from the interaction of 
high energy electrons with the magnetic fields of the 
galaxy, and its signature will be present in both tem¬ 
perature and polarization maps. Giving the depen¬ 
dence of the synchrotron optical depth with frequency, 
the power of synchrotron emission grows with de¬ 
creasing frequency. It is usually modelled using maps 
at 408 MHz [50] and with the WMAP K-band at 23 
GHz [51]. The synchrotron power spectra is well fitted 
using a simple power law for both i and v. The latest 
Planck model [52] is 


Cf = As 




( 22 ) 


where the values of the different parameters are shown 
in Table HI. 


b. Thermal Dust 


1. Instrumental Noise 


The sensitivity of the detectors of a given CMB exper¬ 
iment is finite; thus, a certain noise will be induced in 
the map due to the deconvolution of a Gaussian beam, 
which reads as [45] 






02 ■ 
81n2 


( 21 ) 


Contrarily to synchrotron emission, the power at which 
thermal dust radiates grows with frequency. Planck has 
modelled the dust contamination using a Modified Black 
Body for which To = 19.6 K. The intensity [52] and 
polarization [49] spectra can be written as 



and 


where corresponds to the temperature and polariza¬ 
tion sensitivity of the channel, respectively {X = {T, P}), 
and 9 is the Full Width at Half Maximum (FWHM) of 
the beam. We follow here the specifications for the fu¬ 
ture COrE mission given in Ref. [13], see Table VH of 
Appendix B. 


/^D 

p 


t (y { b,{to) 

2 ^ V^d) \B.i{To) 


2 

(24) 


respectively, where B,^{T) = 2hv^c ^/(e fer—1). The val¬ 
ues of the different parameters are specified in Table HI. 


2. Foregrounds 

Foregrounds, consisting of radio emissions from the 
galaxy and/or other sources at the same frequency 
to that of the CMB signature, will clearly be the 
dominant limiting factors in extracting the cosmological 
information from the maps. In the case of the polarized 
signal, foregrounds are critical as they are orders of 
magnitude higher than the primordial signal in some 
cases. The usual strategy followed to deal with the 
foregrounds is to exploit their spectral dependence. Sev¬ 
eral recent works [46-49] have shown that an accurate 
multifrequency approach to correctly handle foregrounds 
is mandatory. Here we will briefly discuss the physical 
origin of the main foregrounds relevant for the COrE 


3. Statistical Method 


In order to forecast the errors of the different param¬ 
eters we follow the widely used Fisher matrix formal¬ 
ism [53]. The Fisher matrix is defined as the expectation 
value of the second derivative around the maximum of 
the likelihood 


52 £ 

p.. = _(_\|- - 

*2 


(25) 


^ Other two sources of foregrounds are the Anomalous Microwave 
Emission and the Free-Free emission (see Ref. [52] for details 
related to their parametrized power spectra) not discussed here, 
as their impact at the frequency range of interest is negligible. 
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Foreground 

Parameter 

Planck 


As (/tAcmb) 

(4.2 ±0.4) X 10® 


ns (GHz) 

0.408 

Synchrotron 


100 


PSyn 

-3.00 ±0.05 


(^Syn 

-2.5 ±0.02 


Ad (pKqmb) 

40 ±3"- 


UD (GHz) 

353 

Dust 


100 


Pd 

1.51 ±0.01 


an 

-2.4 ±0.02 


Abe (mAcmb) 

247 ± 3 


Abb/Abe 

0.53 ±0.01 


(GHz) 

353 

Dust Polarization 


80 


PI 

1.59 ±0.17 


CtEE 

-0.42 ±0.02 


C^BB 

-0.44 ±0.03 


From Ref. [52], after applying color corrections and conversion 
units. 

From Table 1 of Ref. [49], after applying the color corrections. 

TABLE III. Parameters for the different foregrounds consid¬ 
ered in this study, for both the intensity and polarized emis¬ 
sions. For the intensity signal the models are fitted for £ < 100 
and, for polarization, for 60 < £ < 500. As commonly carried 
out in the literature, we will extrapolate the models to higher 
and lower multipoles for both the intensity and polarization 
spectra. 


where 9i represent a cosmological parameter, and dijid 
represents the fiducial value for the parameter. The 
Cramer-Rao bound ensures that for unbiased estima¬ 
tors the best achievable la error for a given parameter 
marginalized over the other parameters is 

. (26) 

with F~^ the inverse of the Fisher matrix. 


4- Foreground removal 

As argued in the previous section, the main limita¬ 
tion for future CMB observations is the foreground con¬ 
tamination. Among the two polarized foregrounds spec¬ 
ified above, the most dangerous one when measuring the 
tensor-to-scalar ratio r is the galactic dust component, 
as, in general, it gives the largest contribution at the 
Planck and COrE frequencies. 

The issue of foreground removal is a delicate one. 
Many techniques like template cleaning, bayesian esti¬ 
mation, internal linear combination or independent com¬ 
ponent analysis are used for this purpose (see [59] for a 


summary). For example, in Ref. [60], a study forecast¬ 
ing errors on r is performed, without any assumption of 
the properties of the foregrounds. In Refs. [58, 61] the 
errors on the different cosmological parameters are ob¬ 
tained after marginalising over the foregrounds following 
some simple models for their spectra. Here, following 
the approach of [43], we will assume a simple model for 
the foregrounds (see Eqs. (22), (23) and (24)). We shall 
also assume in the following, for simplicity, that the fore¬ 
grounds will be subtracted by a constant amount. Given 
that the Planck mission has achieved a less than 10% 
foreground removal in power, for the COrE mission, due 
to the high number and the high sensitivity of channels 
devoted to the study of the dust, one should expect that 
power could be removed at the 1% level, which is equiv¬ 
alent to set in Eq. (18) (T_f(jz) = 0.01. 

B. Results 

1. Future satellite CMB missions 

In the following, we shall apply the Fisher matrix 
method to the future CMB mission COrE (see Appendix 
A for a consistency check of our method), although simi¬ 
lar results could be obtained for other future CMB satel¬ 
lite experiment. We shall use the 105, 135, 165 and 195 
GHz channels for all the runs, see Table VH in the Ap¬ 
pendix B ®. 

We perform two different analyses. The first one con¬ 
siders no foreground contamination. The second one re¬ 
lies on a 1% foreground subtraction in power {ap = 0.01). 
We assume no delensing on the the B-mode signal. The 
results are shown in Table IV. Comparison between the 
fifth and sixth columns confirm, numerically, the very- 
well known fact that foregrounds will be the major limita¬ 
tion for future CMB missions when extracting the tensor- 
to-scalar-ratio r. 


2. Future constraints on inflationary parameters 

Figures 1, 2 and 3 show, together with the theoreti¬ 
cal predictions and the current constraints from Planck 
measurements, the results of our COrE forecasts for two 
possible values of the tensor-to-scalar ratio (r = 0.1 and 
0.01) and two possible fiducial models. The values of the 
inflationary parameters for the first fiducial model are 
Us = 0.9591, Us = 0.0077 and /3e = 0.0313, which cor¬ 
respond to the best-fit to Planck data. For the second 
fiducial model, which aims to lie within the region cov¬ 
ered by the theoretical models explored here, the values 
are Us = 0.9591, as = —0.0005 and /3s = 0. Tables V 


For recent CMB forecasts see Refs. [54, 55], where however the 
parameters cts and fs where not considered. 











Parameter 

Fiducial 

Planck 2015 

Planck (Fisher forecast) 

COrE ap = O.OU 

COrE (JF = 0 


0.02223 

0.00017 

0.00013 

0.000065 

0.000052 


0.1202 

0.0015 

0.0012 

0.00076 

0.00036 

h 

0.6762 

0.0069 

0.0054 

0.0031 

0.0014 

'^veio 

0.079 

0.019 

0.018 

0.0084 

0.0024 

ln(10i°Ag) 

3.117 

0.037 

0.036 

0.015 

0.0044 

Us 

0.9591 

0.0056 

0.0053 

0.0034 

0.0023 

OLs 

0.0077 

0.011 

0.0077 

0.0040 

0.0036 

Ps 

0.0313 

0.014 

0.019 

0.0088 

0.0065 

r 

0 

< 0.23 

< 0.38 

< 0.016 

< 0.0001 


^ Using TTTEEE + lowP Planck 2015 data. 

^ Using the Planck foreground specifications and the 100, 143 and 217 GHz channels. 

^ Using (jjp = 0.01 and the 105, 135 165 and 195 GHz channels. 

^ Using only resolution noise and the 105, 135 165 and 195 GHz channels. 

TABLE IV. Results for the Fisher Matrix Analysis, 68% CL for all parameters and for r at 95% CL. The fiducial has been 
assumed to be the same for each run as given in the second column. We have used — 2000 and — 500. 


T (Jr 

0.1 0.0096 0.0034 0.0040 0.0086 
0.01 0.0089 0.0034 0.0040 0.0085 


TABLE V. Results from the Fisher matrix analysis. The 
fiducial model corresponds to the Planck 2015 best-fit values 
(us = 0.9591, tts = 0.0077, /3s = 0.0313) and the forecasted 
errors are obtained assuming a foreground removal ap = 0.01 
for the future COrE mission. 


V (Jr ^Tis ^cts 

0.1 0.0096 0.0040 0.0042 0.011 

0.01 0.0088 0.0041 0.0043 0.012 


TABLE VI. As Table V but for a hducial model based on 
the values predicted by the theoretical models explored here 
(ns = 0.9591, as = -0.0005, /3s = 0). 


and VI show the Icr errors on the inflationary parameters 
for these two fiducial models. Notice that the uncertain¬ 
ties on the Us, cts, and j3s barely depend on the fiducial 
value of r, as the tensor-to-scalar ratio is not strongly de¬ 
generate with these parameters. The quantity ar is the 
expected error from the COrE experiment. Notice that 
the error is always larger than the COrE sensitivity limit, 
~ 8 X 10“^ (higher than the target of Ref. [54], which was 
the theoretical milestone of 2 x 10“^). However, the pa¬ 
rameter space, the treatment of the foreground removal 
and the delensing assumptions of future CMB missions 
for the present study and those of Ref. [54] are different. 

From the forecasted errors in Tables V and VI and 
Figure 1 we notice that the measurement of r will not be 
enough to discriminate between the models on the (r, Ug) 
plane. As Figures 2 and 3 show, the forecasted errors on 
the additional parameters considered here, ag and /3g , are 


wider by an order of magnitude or more than the region 
of values for which the most favoured inflationary models 
explored in Section III spread. Thus, there is no hope in 
disentangling between the different models using these 
parameters when the data points to their nominal val¬ 
ues of Ug « —0.0005 and f3g « —10“^. However, notice 
from Figure 3, that if the best-fit value of the /3g parame¬ 
ter arising from future CMB data agrees with its current 
best-fit from Planck measurements, then, this parameter 
could allow to exclude the inflationary models explored 
here at a high confidence level (with the precise signifi¬ 
cance level depending on the particular model under con¬ 
sideration). Finally (see also Ref. [42]), the error bars on 
nt expected from the CorE mission will be ~ 3 — 4 times 
larger than the spread of the slow-roll predicted values 
(as shown in Figure 4) and therefore this parameter does 
not help in disentangling among the possible theoretical 
schemes. 


VI. CONCLUSIONS 

The recent 2015 Planck measurements still allow many 
of the possible theoretical scenarios (as quadratic-like, 
Higgs-like and Hilltop models) as the underlying infla¬ 
tionary mechanism. A firm confirmation of the inflation¬ 
ary paradigm would require a detection of the primordial 
gravitational wave signal. However, in order to single 
out a theoretical model, the usual two slow-roll parame¬ 
ters, that is, the scalar spectral index Ug and the tensor- 
to-scalar ratio r, may not be sufficient. The reason is 
due to the fact that the (rig, r) plane appears to be un¬ 
evenly filled, with a potentially forbidden zone and other 
highly populated regions in which mostly all the theo¬ 
retical predictions lie. In this regard, we have explored 
the discriminating power of two other observables, the 
running ag and the running thereof /3g. Our analyses 
of Planck temperature and polarization data show that 
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FIG. 1. 68% and 95% CL allowed contours from Planck TT and lowP and from Planck TTTEEE plus lowP in the (ns,r) 
plane. We show as well the forecasted 68% and 95% CL contours from the CMB mission COrE considering a 1% foreground 
removal in power (ctf = 0.01) and perfect foreground subtraction {(jf = 0). Notice that if the residual foreground emission 
is only removed to 1% level, the future COrE mission may no disentangle among the different models in the (ns,r) plane. 
However, the level at which the foregrounds could be removed may be lower than 1%, but this will only be known once future 
measurements are performed. 
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FIG. 2. As Figure 1 but in the (tIs^ols) plane, assuming a 1% foreground subtraction in power. 
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FIG. 3. As Figure 1 but in the {ns,Ps) plane, assuming a 1% foreground subtraction in power. If nature has chosen a value of 
P close to the current best-fit value from Planck, then the inflationary models considered here could be excluded by at a high 
confidence level by future CMB observations. 


the current errors on the former two quantities are large, 
and therefore they do not help in discarding some of the 
possibilities, even if the present mean value of Ps lies 
2a above its predictions in the most favoured inflation¬ 
ary models explored here. However, future CMB mea¬ 
surements, such as the COrE mission, have the poten¬ 
tial to rule-out some theoretical possibilities at a much 
higher significance, provided the best-fit values for these 
additional parameters do not change significantly from 
their current estimates. Our forecasts (which rely on 
both a simple model for foregrounds and assume an ad- 
hoc 1% foreground removal) show that COrE may help 
enormously in unraveling the inflationary mechanism via 
its measurement of /3s, especially if Nature has chosen a 
value of r > 0.005, which is close to the sensitivity limit 
found in this study. Other complementary information 
concerning /3s and/or as, as those coming from future 
planned galaxy surveys [56] (for instance, the SPHEREX 
project [57]), could significantly improve the sensitivities 
forecasted here. 
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VII. APPENDIX 

A. Consistency of the Fisher Method 

We test the validity of our method by computing our 
Eisher matrix forecast for the complete Planck mission 
and comparing our results to those obtained by Planck 
measurements. For that purpose, we shall use the 100, 
143 and 217 GHz channels of Planck with its accounted 
foreground removal as shown in Table HI, and following 
the specifications detailed in Table VH of the Appendix 
B. From the results depicted in Table IV, notice that 
there is an excellent agreement between the forecasted 
parameter errors and the errors quoted by the Planck 
collaboration, with the differences always below the 20% 
level. In addition, we have verified that the correlations 
between the cosmological parameters are well accounted 
for. 
























11 


10 " 


q 

"cS 

_ro 

ro 

o 

c/3 

I 

o 

I 

o 

CO 


10 - 


10 - 



Natural inflation 


Higgs-like model 


Hiltop quadratic model 

— 

Nonminimal chaotic model 

— 

n + n^- 

nflation 

— 

V cc (p‘^ 


— 

V <xp 


• 

AT = 50 


• 

V = 60 



-0.040 -0.035 -0.030 -0.025 -0.020 -0.015 

Tensor spectral index (n^ 


- 0.010 


-0.005 


0.000 


FIG. 4. Theoretical predictions in the (nt, r) plane for the most favoured inflationary models studied here. We do not illustrate 
the current nor the forecasted contours, as the size of the region depicted here lies well within their Ict range, showing that nt 
can not help much in distinguishing among these theoretical possibilities (see the recent [42]). 


B. CMB Mission specifications 

Table VII shows the values used for the sensitivity of 
Planck and COrE missions, as in Ref. [13]. 
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